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ABSTRACT 

We present the results of the kinematic study of the extended gas in a sample of 10 
high redshift radio galaxies {z ~2.5) based on high signal to noise Keck II and VLT 
long slit spectroscopy. In addition to the typical high surface brightness kinematically 
perturbed regions (FWHM and velocity shifts >1000 km s~^), we find in all objects 
giant low surface brightness halos which show quieter kinematics with typical emission 
line FWHM and velocity shifts of ^several hundred km s~^. 

The giant halos often extend for more than 100 kpc and sometimes beyond the 
radio structures. They emit lines other than Lya (CIV, Hell and NV in some cases), 
typically found in the spectra of high redshift active galaxies. Continuum is also often 
detected. The halos are enriched with heavy elements at tens of kpc from the active 
nucleus. Typical Lya luminosities and surface brightness (within the slit) are in the 
range lO'''^"'*'* erg s~^ and several x 10~^^ *° erg cm~^ s~^ arcsec"^ respectively. 
Estimated densities are in the range ^^17-150 cm^'^. The quasar continuum is the 
dominant source of ionization of the quiescent halos along the radio axis. The implied 
total quasar ionizing luminosities are in the range ^several x lO'^^-lO'^'^ erg s ^, in 
the same range as radio loud quasars at comparable redshift. 

The detection of giant quiescent halos in all objects suggests that they could 
be a common ingredient of high redshift radio galaxies. The radio galaxies seem to 
be embedded within the halos. The nature, origin of the halos and the cosmological 
implications are also discussed. 
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1 INTRODUCTION 

Our current belief is that the hosts of powerful radio sources 
in the distant Universe are destined to become the giant 
ellipticals of today, the most massive galaxies we know (e.g. 
Lilly & Longair 1984, McLure et al. 1999). They inhabit rich 
environments that are believed to be galaxy clusters in the 
process of formation (e.g. Pentericci et al. 2000, Venemans 
et al. 2002). While some distant powerful radio sources may 
have commenced their formation at very high redshift, the 
process of assembly is very active at z ~2-3. This corres- 
ponds to the epoch when luminous quasars appear to have 
had their maximum space density and the star formation 
rate of the Universe was two orders of magnitude higher 
than at present (e.g. Pei 1995, Shaver et al. 1996, Fan 2001). 

Our group has recently proved that solar or superso- 
lar metallicities are common in the extended gas (tens of 
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kpc) of radio galaxies at z ~2.5 (Vernet et al. 2001, Villar- 
Martm et al. 2001). The relatively high metallicities imply 
that star formation in these objects has been rapid and in- 
tense, enabling a fast chemical enrichment of the gas. Since 
stellar population analysis of giant ellipticals at low redshift 
indicates that most stars have been formed during a brief 
epoch with high star formation rate, this further supports 
the idea that powerful high redshift radio galaxies are pro- 
genitors of the giant ellipticals of today. Star formation is 
also suggested by the detection at submm wavelengths of 
some high redshift radio galaxies (e.g. Archibald et al. 2001, 
Reuland et al. 2003a). Although the nature of the heating 
mechanism is controversial, some authors have argued that 
the dust is heated by starbursts, rather than the hidden 
quasar nuclei (e.g. Tadhunter et al. 2002). By looking at 
high redshift radio galaxies [z > 2-3, hereafter HzRG) we 
are therefore likely to be witnessing the formation process 
of giant elliptical galaxies. 

HzRG are surrounded by giant nebulae which often ex- 
tend for more than one hundred kpc (e.g. McCarthy et al. 
1990a, Reuland et al. 2003b) and sometimes beyond the 
radio structures (e.g. Eales et al. 1993, Kurk et al. 2001, 
Maxfield et al. 2002). One of the great difficulties in the in- 
terpretation of the optical and near infrared data of HzRG 
is due to the fact that the observed morphological, kine- 
matic and ionization properties are strongly influenced by 
the nuclear activity and the extreme phenomena associated 
with the relativistic jets. Most morphological and kinematic 
studies of HzRG have been focussed on the high surface 
brightness regions. These are clumpy, irregular (with fea- 
tures such as filaments, plumes, ionization cones, e.g. Reu- 
land et al. 2003b) and often aligned with the radio axis 
(Chambers, Miley & van Breugel 1987, McCarthy et al. 
1987). They are characterized by extreme kinematics, with 
measured FWHM and velocity shifts >1000 km s~^ (e.g. 
McCarthy, Baum & Spinrad 1996, Villar-Martm, Binette & 
Fosbury 1999) , compared to values of a few hundred in low 
redshift radio galaxies (e.g. Tadhunter, Fosbury & Quinn 
1989, Baum, Heckman & van Breugel 1990). It is likely 
that a perturbing mechanism such as shocks induced during 
jet/gas interactions ( |van Ojik et al. 1997^ is responsible for 
such extreme kinematics. 

However, there is evidence in two HzRG for giant (dis- 
tance to the nucleus > 50 kpc^) low surface brightness ha- 
los (LSBHs hereafter) of ionized gas with quiet kinematics: 
1243-1-036 at z =3.6 (van Ojik et al. 1996, v096 hereafter) 
and 0828-1-193 at z =2.6 (Villar-Martm et al. 2002, VM02 
hereafter) . 

These quiescent halos might be residual gas from which 
the galaxy started to form (VM02). They offer opportunities 
to study the ionizing photon luminosity of the quasar and 
with sufficient line measurements the chemical composition 
of the halo gas. If the gas kinematics is gravitational in origin 
the halos can be used as mass tracers of radio galaxies at 
different redshifts. They would provide critical information 
to constrain models of galaxy formation and evolution. Also, 
the metallicities of the giant halos will set tight constraints 
on their chemical enrichment history and, ultimately, on the 
formation process of the galaxy. 



1 We assume nA=0.7, nAf=0.3, Ho=65 km s"! Mpc"! 



We present in this paper the results of a research project 
whose goal is to search for kinematically unperturbed halos 
in HzRG and study their kinematic and ionization proper- 
ties. An important advantage of the current set of data is 
that we can use Hell (and CIV in some cases) to constrain 
the kinematic properties of the gas, rather than using only 
Lya, which is very sensitive to dust/gas absorption. 

The paper is organized as follows: The observations and 
data analysis are described in §2. General results and results 
on individual objects are presented in §3. The ionizing mech- 
anism, kinematic patterns, nature of the giant quiescent ha- 
los and cosmological implications are discussed in §4. The 
summary and conclusions are presented in §5. 



2 OBSERVATIONS AND DATA ANALYSIS 
2.1 The data 

All the sources were selected from the ultra-steep spectrum 
(USS) radio galaxy survey (see e.g. Rottgering et al. 1994) 
with redshift z >2.2. 

The spectra for all the objects (see Table 1) except 
2104-242 were obtained with the Low Resolution Imaging 
Spectrometer (LRIS, Oke et al. 1995) with its polarimeter 
(Goodrich, Cohen, Putney 1995) at the Keck II 10 m teles- 
cope from July 1997 till July 1999 under subarcsec seeing 
conditions (seeing varying between ~0.7 and 1.1 arcsec). The 
observing log is presented in Vernet et al. (2001). 

The LRIS detector is a TEK 2048^ CCD with 24 ^.m 
pixels giving a scale of 0.214 arcsec pixel"^. We used a 
300 line mm~^ grating which provides a dispersion of 2.4 
A pixel"^ and a spectral range Xobs ~3900-9000 A. The 
slit width was always 1 arcsec, except for 4C-00.54 and 
4C-I-23.56 for which we used a 1.5 arcsec slit. The slit was 
always oriented along the radio axis given by Rottgering et 
al. 1994 and Carilh et al. 1997 (see Table 1). 

Each individual frame was calibrated in wavelength and 
corrected for slit curvature using arc spectra. This initial 
wavelength calibration was refined using strong sky lines. 
The FWHM of the instrumental profile (IP, Table 1) was 
determined from unblended sky lines. These values corres- 
pond to a range in resolution of 320-610 km s~^ at the 
redshifted HeIIA1640 wavelength. Variations from object to 
object are probably due to mechanical instabilities within 
the instrument causing variation on the line spread function. 

In order to compensate for the effects of fiexure in 
the instrument, we extracted from every frame a section 
around Lya and NVA1240 and one around CIVA1550 and 
HeIIA1640. The individual sub-frames corresponding to the 
same spectral range were spatially aligned (using the con- 
tinuum centroid) and combined. The two final frames (one 
for Lya and NV and another one for CIV and Hell) were 
also spatially aligned, assuming that the continuum centroid 
has the same spatial position at all wavelengths. This is rea- 
sonable in this small spectral range, taking also into account 
that the spatial resolution element has a large physical size 
(~ 8-9 kpc). For a more detailed description of the observa- 
tions and data reduction see Vernet et al. (2001). 

The observations for 2104-242 were carried out in ser- 
vice mode on 1999 September 2-5 with FORSl on the 8.2 m 
VLT Antu telescope (ESO-Chile). The 600B grism was used 
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with a 1 arcsec wide slit. The slit was positioned along the 
brightest components and the filamentary structure at a po- 
sition angle of 2° north through east. The exposure time was 
10800 sec. The seeing during the observations was 1 arcsec 
and conditions were photometric. A 2 x 2 readout binning 
was used in order to increase the signal-to-noise ratio (S/N). 
The resultant spectral resolution was 6 A (FWHM) and the 
spatial scale was 0.4 arcsec pixel"^. For a detailed descrip- 
tion of the data reduction see Overzier et al. (2001). 

Radio maps, optical and near infrared HST archive ima- 
ges were also used and registered with the Keck 2D spectra 
for comparison between the structures observed at different 
wavelengths. The radio observations and data are described 
by CariUi et al. (1997) and the HST WFPC2 and NICMOS 
images are described in Pentericci et al. (1999 and 2001) 
respectively. 

2.2 Registration 

Several methods of registration between the radio maps and 
the 2D spectra were used, depending on the availability of 
HST images and on the presence of a radio core. 

For objects presenting a radio core and for which 
we have both WFPC2 and NICMOS images (0211-122, 
0828+193, 1410-001 and 2104-242), the IR peak was assu- 
med to coincide with the radio core. The WFPC2 image was 
then registered with the NICMOS image using field stars. 
Finally the 2D spectrum continuum centroid was aligned 
with the WFPC2 image peak. The total uncertainty with 
this method is ~0.3 arcsec. The error might be larger: since 
the AGN is expected to be obscured, the near infrared peak 
probably does not mark its position and therefore, it might 
be spatially shifted relative to the radio core. 

2105-1-236 is the only object presenting a radio core for 
which we only have a WFPC2 image. The obvious thing to 
do would have been to align the radio core with the peak of 
the WFPC2 image. However, the available astrometry places 
the radio core on top of a secondary peak of the WFPC2 
image, close to the Lya emission vertex in Chambers et al. 
H1996|l ground-based narrow band image. We used this re- 
gistration point instead. 

Three objects (0731+438, 1558-003 and 1931+480) 
have a radio core but HST images are not available. We sim- 
ply aligned the radio core with the peak of the continuum 
in our 2D spectra. The total uncertainty with this method 
is ~0.2 arcsec. As above, this error is likely to be larger 
since the continuum centroid might not mark the position 
of the active nucleus. (Notice that for 1931+480 we have 
used a different radio core identification than Chambers et 
al. 1996). 

Two objects do not have a radio core. For 1809+407 
the radio map and optical spectrum were aligned using the 
astrometry from Chambers et al. 1199611 . We assumed that 
the optical continuum centroid in the Keck spectrum coin- 
cides with the optical identification by Chambers et al. The 
uncertainty of this method is ~0.5 arcsec. For 0943-242, we 
determined the astrometry of field stars in the WFPC2 ima- 
ge using the GSC2.2 catalogue (accuracy of ~0.2-0.3 arcsec) 
and then aligned it with the radio map. The continuum cen- 
troid in the 2D spectrum was then aligned with the peak of 
the image. The uncertainty of this method is ~0.4-0.5 arc- 
sec. 



2.3 Data analysis method 

The 2D spectra were divided into several apertures along 
the spatial direction and a spectrum was extracted for each 
one. The apertures were selected so that the gas has simi- 
lar kinematic properties across the spatial extent of a given 
aperture and to obtain enough S/N ratio in Lya, CIV and 
Hell (if possible) to fit the line profiles. Separate apertures 
were defined for the gas beyond the radio structures when 
detected. 

The Lya and Hell spectral profiles were analysed and 
fitted in each aperture with one or more Gaussians depen- 
ding on the quality of the fit. For each kinematic compo- 
nent the FWHM (corrected for instrumental broadening in 
quadrature) and the velocity shift relative to the Hell line 
at the continuum centroid were calculated. Hell has been 
the most important line in our study because it is a strong 
single (rather than a doublet) non resonant line. Lya, in 
spite of its strong sensitivity to absorption by gas and dust 
has also often provided very valuable kinematic information. 
CIV is difficult to constrain due to the presence of the dou- 
blet components. Results for this line are only shown when 
they provide convincing evidence for the presence of low ve- 
locity gas, a signature of the quiescent halos. 

Sometimes different fits produced results of similar 
quality. In most cases it was possible to discriminate by 
comparing with the results from (a) adjacent apertures (b) 
another emission line detected in the same aperture. As an 
example, it was sometimes found that for a given spectrum 
one Gaussian component produced as good a fit to the Hell 
profile as two Gaussians components. The profile of Lya in 
the same aperture and the fits to the Hell lines in the ad- 
jacent apertures provided further constraints that sugges- 
ted the use of either one or two components. When such 
discrimination was difficult, no spectral decomposition was 
attempted. 

Any fit which produced lines with observed FWHMots 
such that 

FWHMobs < IP - |AIP| - |AFWHMo6s 

was rejected, where AIP is the uncertainty on the 
FWHM of the IP and AFWHM^b, is the error on FWHM^i,. 
A kinematic component was considered to be unresolved 
when FWHMoijs was in the range 

IP - |AIP| - |AFWHMo6,| < FWHMob, < IP + |AIP| 
+ |AFWHMoi„|. 

In such cases an upper limit for the intrinsic FWHM 
was set by estimating the intrinsic FWHMi„t for which 
FWHMoSs = IP + |AIP| + |AFWHMo6s|. 

The formal uncertainty in the central wavelength of an 
emission line is a combination of systematic uncertainties 
in the wavelength calibration (estimated to be <1 A), the 
signal to noise ratio of the line and the errors on the Gaus- 
sian fits. The formal uncertainty in the deconvolved width 
(FWHMint) is a combination of uncertainties in the resolu- 
tion, the signal to noise ratio of the line and the errors on 
the multiple Gaussian fits. 

The real errors are likely to be larger because the seeing 
disk might have been smaller than the slit width. Although 
we are confident that the seeing was always 0.7 arcsec, 
the assumption that the objects filled the slit introduces 
additional uncertainties which will be discussed in §3.2. Our 
conclusions are not seriously affected by this assumption. 
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lAU 


4C 


Redshift 


Exp. 
sec 


PA 


IP 

A 


HST 
images 


1809+407 


4C+40.36 


2.265 


10742 


82 


9.8+1.2 


No 


0211-122 




2.340 


28580 


104 


9.5+1.2 


WFPC2/NICMOS 


1931+480 


4C+48.48 


2.343 


12000 


50 


11.2+1.2 


No 


1410-001 


4C-00.54 


2.360 


16240 


134 


14.3+1.2 


WFPC2/NICMOS 


0731+438 


4C+43.15 


2.429 


22720 


12 


10.5+1.2 


No 


2105+236 


4C+23.56 


2.479 


22266 


47 


13.6+1.2 


WFPC2 


2104-242 




2.491 


10800 


2 


6.0+0.5 


WFPC2/NICMOS 


1558-003 




2.527 


11400 


72 


11.0+1.2 


No 


0828+193 




2.572 


18000 


44 


10.5+1.2 


WFPC2/NICMOS 


0943-242 




2.922 


13800 


73 


10.8+1.2 


WFPC2 



Table 1. The radio galaxy sample. Objects have been organized in order of increasing redshift. The exposure times, slit position angles 
(PA) and instrumental profiles FWHM (IP) are also shown. Results for 0828+193 have been presented in Villar-Martfn et al. 2002. The 
last column shows the HST images used in the regitration process between the radio/HST images and Keck spectra (§2.2) 



Name 


Rrad 


Rraax 


FWHM 


Vs 






Detected 












10-17 


10« 


Lines 




kpc 


kpc 


km s^-' 


km s^^ 


erg s"! cm" ^ arcsec— 2 


erg s^i 




1809+407 


27 


40 


500-700 


400 


2.0 


1.6 


Lya, NV (?), CIV, Hell 


0211-122 


62 


61 


<400-700 


200 


1.6 


1.4 


Lya, NV, CIV, Hell 


1931+480 


110 


53 


450-720 


450 


7.1 


24.7 


Lya, CIV, Hell 


1410-001 


113 


62 


<472-800 


500 


2.2 


10.5 


Lya, CIV, Hell 


0731+438 


42 


35 


600 


<100 


48.9 


11.6 


Lya, CIV, Hell 


2105+236 


170 


61 


500-650 


570 


13.6 


7.4 


Lya, Hell 


2104^242 


135 


65 


<265 


250 


2.5 


0.9 


Lya, HeII(?) 


1558-003 


61 


69 


550-850 


550 


2.3 


13.4 


Lya, CIV, Hell 


0828+193 


64 


69 


<300 


600 


2.4 


10.2 


Lya, CIV, Hell 


0943-242 


17 


59 


400-600 


450 


1.3 


1.2 


Lya, NV, CIV, Hell 


1243+036 


31 


78 


250 


450 




3.0 


Lya (at least) 



Table 2. Main properties of the quiescent giant halos. Rmax is the maximum extension measured from the position of the continuum 
centroid using the Lya line. Rrad is the extension of the radio source measured from the position of the continuum centroid on the same 
side as Rmax- Column (4) gives the range of maximum and minimum FWHM values measured across the quiescent halos. Vg is the 
maximum velocity shift measured across the halo. Both FWHM and Vs refer to the Hell line except for 2104^242 for which Lya was 
used (see text §3). Surface brightness (SB) and luminosities (within the slit) refer to Lya, except for 0943-242 for which Hell was used 
(see text §3). The results for 1243+036 and 0828+193 were taken from v096 and VM02 respectively. 



3 RESULTS 

Figures 1 to 9 show the data and results of the kinematic 
study. Each figure corresponds to one radio galaxy and con- 
sists of 3 panels: the top panel shows the overlay of the ra- 
dio map and the WFPC2 image (when available) spatially 
aligned with the 2D Lya, NV, CIV and Hell spectra. The 
apertures selected for the kinematic analysis are indicated. 
The bottom left panel shows the ID spectra extracted from 
the selected apertures and the bottom right panel shows the 
spatial variation of the kinematic properties (velocity shift 
and FWHM) of the difi^erent components revealed by the 
spectral fits. When two kinematic components were isolated, 
these are distinguished with open (narrow components) and 
solid symbols (broad components). Circles are used for Lya, 
triangles for Hell and squares for CIV. The dashed vertical 
lines indicate the maximum extent of the radio structures. 
The main results of the kinematics analysis are: 
1) All 10 HzRG (including 0828+193, VM02) are 
associated with giant halos of quiescent gas with 
FWHM(HeII)<850 km s"^ and velocity shifts across the 
nebulae <600 km . 



2) In addition, 8 out of 10 objects contain kinemati- 
cally perturbed gas. 7 of these objects show a broad compo- 
nent with FWHM>1500 km s"^ and high surface brightness. 
2104-242 (which was observed at higher resolution) shows 
signs of kinematic perturbance in the form of split narrow 
components. The perturbed gas is usually located inside the 
radio structures. 

We show in Table 2 some properties of the quiescent 
LSBHs. We have included 0828+193 (VM02) and 1243+036 
(v096) for comparison. Rmax is the maximum extension 
measured from the position of the continuum centroid us- 
ing Lya (the brightest line). Halo sizes axe in the range 
Rmax ~35-80 kpc^ and they extend sometimes beyond the 
radio structures (it is possible that radio maps with greater 
sensitivity might show that the radio plasma extends fur- 
ther out). The giant halos emit lines other than Lya (CIV, 
Hell and NV in some cases), typically found in high red- 
shift active galaxies. The emission line spectra are typical of 

^ In the Cosmology we have adopted, Rmax ~78 kpc for the halo 
detected by v096 in 1243+036. 
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active galaxies. Continuum is also often detected. In most 
cases the quiescent gas is found not only in the outer parts 
of the objects but across the high surface brightness regions 
as well. 

The kinematic parameters FWHM (maximum and mi- 
nimum values across the halo) and Vs (maximum velocity 
shift measured across the halo) in Table 2 refer to Hell, 
which is a more reliable kinematic tracer than Lya. Lya was 
used only for 2142-242, since this was the only line detected 
from the quiescent halo (§3.1.7). 

The Lya emission of the quiescent halos is clearly de- 
tected in the outer regions of all objects except 0731-1-193 
and 2105-1-236 (see top panels Fig. 1-9). The surface bright- 
ness values presented in Table 2 were measured using these 
regions. It is important to note that the surface brightness 
varies spatially and it is often higher in the inner apertures. 
For 0731+193 and 2105-1-236 the surface brightness was esti- 
mated using the flux of the narrow component isolated with 
the kinematic analysis and the area of the apertures where 
this component was detected (see§3.1.5). 

The Lya luminosities were calculated integrating the 
line flux from the quiescent gas across all the apertures 
where this was isolated. For 0943-242 it was not possible 
to isolate a narrow component using Lya on the high sur- 
face brightness regions (see §3.1). In this case we show the 
luminosity of the Hell narrow line. 

3.1 Notes on individual objects 

3.1.1 1809+407 

Narrow band Ha+[NII] imaging obtained by Egami et al. 
H2003|l using NICMOS on HST shows that the emission hne 
nebulosity is misaligned from the radio axis. Our 2D Lya 
Keck spectrum (Fig. 1) shows that in addition to the high 
surface brightness regions, a LSBH with apparently quieter 
kinematics (narrower Lya) is detected extending ~9 arcsec 
(~80 kpc) and beyond the radio structures. 

Spectra were extracted from 5 apertures (Fig. 1). The 
fits to Hell reveal the presence of a narrow component with 
FWHM in the range ~500-700 km s~^ and maximum ve- 
locity shift across the nebula of ~400 km s^^. These results 
are consistent within the errors with those obtained using 
Lya. Our results are in very good agreement with Egami 
et al. ()2003|l . These authors found that the integrated near 
infrared (optical rest frame) emission lines can be characteri- 
zed by two components of FWHM~560 km s~^ and 1670 km 
s~^ respectively and shifted by ~-470 km s~^. 

3.1.2 0211-122 

The 2D line spectra (Fig. 2, top) show evidence for a low 
surface brightness halo which extends in Lya for ~13 arcsec 
or ~114 kpc. The apparent kinematics of the LSBH seems 
rather uniform across its entire extension and quieter than 
in the high surface brightness regions. 

Spectra were extracted from six apertures (Fig. 2). Both 
CIV and Hell (except aperture 3, see below) are very well 
fitted with single Gaussian profiles. This is not always the 
case for Lya, which is heavily absorbed in some positions 
l |van Ojik et al. 1994| l and therefore not a reliable kinematic 
tracer in this object. The fits to the line profiles confirm that 



the kinematics is quieter in the outer parts of the object. 
Both Hell and CIV are unresolved (FWHM<400 km s"^) 
and the velocity shift (< 100 km s~^) between these two 
lines is small. 

Lya is broad (FWHM~ 1000- 1200 km s~^) in the inner 
apertures (ap. 3 and 5). . Similar highly perturbed kine- 
matics is not apparent in Hell, which shows FWHM~700 
km s~^. The CIV (AA1548,1551) profile is also consistent 
with two components of FWHM~700 km s~^ and this is 
also the case for the NV (AA 1239,1243) doublet (ap. 4) 
when broad underlying wings are taken into account. We 
therefore conclude that the gas in the inner apertures emits 
lines with FWHM~700 km s"^, except for Lya. The maxi- 
mum velocity shift for the Hell line across the entire object 
is ~200 km s"\ 

It is interesting to note that NV is very strong relative 
to Hell and CIV at a large distance (~45 kpc) from the 
nuclear region (see spectrum ap. 6, Fig. 2). This will have 
implications on the chemical abundances (see §4). 

3.1.3 1931+480 

The Lya image by Chambers et al. (1996) shows a centrally 
peaked nebula which extends for ~10 arcsec aligned with 
the radio axis. The Keck spectrum (Fig. 3) shows that the 
Lya emission extends for ~13.5 arcsec (~118 kpc). The low 
surface brightness emission shows apparently narrower Lya 
in the outer regions of the object. Spectra were extracted 
from five apertures. 

The fits to Hell reveal a narrow component in all 
apertures (2 to 5) with FWHM in the range ~450-720 
km s~^ and maximum Vs ~450 km s~^ across the ob- 
ject. The K-l-H image of these object shows four distinct 
clumps llCarson et al. 200111 . Near infrared Keckll spec- 
troscopy of the individual clumps reveals that the FWHM 
of [OIII]A5007 is ~525, 390, 840 and 570 km s"^ for the in- 
dividual clumps. These results are in good agreement with 
ours. 

3.1.4 1410-001 

There is a strong misalignment on this object between the 
radio axis and the optical structures as seen in the HST 
images (~45 deg, Pentericci et al. 2001). The slit was located 
along the radio axis, therefore the regions within the slit lie 
far beyond the HST structures (except maybe in the inner 
2-3 arcsec). 

The 2D Lya spectrum shows low surface brightness 
emission in the outer parts of the object with apparently 
quieter kinematics (Fig. 4). 

Spectra were extracted from seven apertures. Lya (ap. 
7), Hell (ap. 5 and 6) and CIV (ap. 6) are very na- 
rrow (unresolved) in the outer parts of the object with 
FWHM(HeII)<472 km s'^. Hell is broader (~700-800 km 
s~^) in the inner apertures. The maximum velocity shift 
across the object is ~500 km s~^. A narrower component is 
suggested by the very narrow peak found on top of the CIV 
and Hell profiles at some spatial positions. 

^ Van Ojik measured FWHM(Lya) ~ 1700 km in the central 
2 arcsec. The authors worked at 24 A spectral resolution. 
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3.1.5 0731+438 

Narrow band Ha+[NII] images obtained with the Subaru 
telescope (Motohara et al. 2000) show diffuse line emission 
which extends out to -^3.3 arcsec from the center. The mor- 
phology is aligned with the axis defined by the radio hot 
spots. The authors interpret the observed morphology as 
evidence for the ionization cones in this radio galaxy. The 
Keck spectrum (Fig. 5) shows that Lya extends for ~11 
arcsec (^^96 kpc). 

Three apertures were selected for the kinematic study 
(Fig. 5). A narrow component with FWHM~600 km s~^ and 
Vs <100 km s~^ is detected in apertures 2 (Lya and Hell) 
and 3 (Lya, CIV and Hell). The emission from aperture 1 
corresponds to a separate component clearly seen in the 2D 
spectra which emits broad lines (FWHM~1000 km s~^). We 
have not considered this component as part of the quiescent 
LSBHs. The results obtained with the three emission lines 
are in good agreement. 

3.1.6 2105+236 

Lya narrow band imaging shows two opposing ionization 
cones with ~90° opening angles (Knopp & Chambers 1997) . 
This is one of the rare cases among high redshift radio 
galaxies with clear ionization cones. The 2D Keck spectrum 
(Fig. 6, top) shows emission from these two spatial com- 
ponents in all emission lines. Contrary to the rest of the 
objects in this sample, the 2D spectrum does not show clear 
evidence for a quiescent LSBH. 

Spectra were extracted from six apertures (Fig 6) . The 
fits both to Lya and Hell are rather complex in this object 
due to the presence of an underlying very broad component 
(probably scattered light, Vernet et al. 2001) detected across 
the whole extension of the SW component, which makes the 
setting of the continuum level uncertain. We have attempted 
to fit the main peak of the emission line as accurately as 
possible, rather than the broad underlying wings, which are 
difficult to reproduce. The errorbars in Fig. 6 (bottom right 
panels) account for the different possibilities considered as 
the continuum level for the fits. This uncertainty is large for 
the high velocity component, but it does not affect seriously 
the narrow component, at least in HeH. 

HeH is narrow in all apertures (2 to 6) with 
FWHM~500-650 km s"^ and the maximum velocity shift 
across the nebula is ~570 km s~^. CIV is also narrow in ap. 
5 with FWHM~600 km s~^ and similar velocity shift. 

The fit to LyQ reveals the presence of a broad com- 
ponent (FWHM~1800 km s"^) and a narrower component 
with FWHM~500-1100 km s~^ Except in ap. 5, both the 
FWHM and velocity shift of this component are consistent 
within the errors with the Hell measurements. 

3.1.7 2104-242 

The 2D VLT spectra (Fig. 7) show emission from two spa- 
tial components that whore described for the first time by 
McCarthy et al. |1990b ). These two clumps are between and 
aUgned with the radio lobes. The Lya emission extends for 
~15 arcsec (~130 kpc), well beyond the structures seen in 
the Lya image (Pentericci et al. 1999). Very narrow Lya is 
detected towards the north-east. 



Spectra were extracted from 6 apertures. The kinematic 
structure (see also Overzier et al. 2001) of this object is 
very complex, but recall that the spectrum has higher spec- 
tral resolution than for the rest of the sample. Two narrow 
(FWHM(HeII)<200 and ~400 km s"^ respectively) compo- 
nents are found in the southern clump, which is associated 
with the narrow filament more than 2 arcsec long detected 
in HST images (Pentericci et al. 1999). 

An interesting feature is the very narrow Lya emission 
(FWHM<265 km s~^) detected in the outer regions of the 
object (ap. 1 and 2). The fit to the Lya profile in ap. 6 reveals 
a similar component (similar FWHM and velocity shift) also 
in the outer regions of the object to the south. Absorption 
of Lya photons could be responsible for the narrowness of 
the line, but the similarity of the kinematic structure of all 
the lines shown by the 2D spectra and the large strength of 
Lya relative to other lines indicates that absorption is not 
important in this object. 

It is uncertain whether the outer very narrow line LSBH 
is also detected across the high surface brightness regions. 
Our spectra do not reveal a component of similar FWHM 
and velocity shift but it might be concealed in the two com- 
ponents of the southern clump. Higher resolution (2.8 A) 
spectroscopy reveals three components (rather than 2) with 
velocity shifts of 140, 670 and 1270 km s~^ respectively (Pen- 
tericci et al. 2001; the authors do not provide FWHM va- 
lues). One of these three components might correspond to 
the outer LSBH. In Table 2, we have assumed that the quies- 
cent LSBH is the very narrow line emitting gas in the outer 
parts of the object. 



3.1.8 1558-003 

Low surface brightness Lya emission is detected beyond the 
radio structures extending for ~14 arcsec (~120 kpc) with 
apparently narrower lines (Fig. 8, top). Spectra were ex- 
tracted from 6 apertures. The fits reveal a narrow compo- 
nent with FWHM in the range <550-630 km s"^ in Lya (ap. 
1), Hell (ap.2) and CIV (ap. 5). The maximum velocity shift 
for Hell across the nebula is Vs ~550 km s~^. 

Hell is broader in the central regions (FWHM~750-850 
km s~^) but the presence of a narrow component is sugges- 
ted by the very narrow peak found on top of the CIV and 
Hell profiles (ap. 4). Although the errors are large (~40% for 
the FWHMobs of the narrow component) the spectral profile 
of Hell in ap. 3 and 4 is well reproduced by a narrow unre- 
solved component with FWHM<410 km s"^ and K ~-|-230 
and -1-340 km s~^ respectively, and a broad component of 
FWHM ~1000 km s"\ (Because of the large errors, we have 
ignored the results of these fits in Table 2). 

It is interesting to note that broad Lya (FWHM~1200 
km s~^) is found beyond the radio structures. This was also 
reported by Villar-Martm, Binette & Fosbury (1999). 

3.1.9 0943-242 

Clear evidence for a low surface brightness very extended 
halo with apparently quiet kinematics is found in the 2D 
spectrum of all (Lya, NV, CIV and Hell) lines (Fig. 9, top). 
This gas extends for ~8 arcsec (67 kpc) from the continuum 
centroid, well beyond the radio structures. 
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Spectra were extracted from 5 apertures. The double 
peaked profile found in some of them is due to absorption 
(e.g. Rottgering et al. 1995). Since the effect is so dramatic 
on the profile of the line, we have focussed the kinematic 
analysis on the Hell line instead. 

A narrow component of FWHM~400-600 km s~^ is 
found in Hell in all apertures with maximum velocity shift 
~450 km s~^ across the nebula. The results for CIV are in 
very good agreement. An interesting feature is the strength 
of the NV line in the quiescent halo relative to other lines 
(Fig. 10). 



3.2 Seeing disk and slit width: Uncertainties 

So far we have assumed that the objects filled the slit. In 
some cases, the seeing disk was smaller than the slit width 
and our assumption introduces additional uncertainties. On 
one hand, if the objects are clumpy, the image of a spatially 
unresolved clump may be smaller than the slit width and the 
instrumental profile at that particular position is narrower 
than the profile measured with the sky or arc lines. This 
would lead to an underestimation of the intrinsic FWHM 
of the emission lines. On the other hand, it is possible that 
some clumps may be displaced in the dispersion direction 
relative to the center of the slit. The net effect will be errors 
in the velocity fields since these may partially reflect the 
distribution of the emitting gas within the slit. 

We do not expect these uncertainties to affect the re- 
sults on 0731+438, 0828+193, 1809+407, 1931+480 and 
2104-242. For 0731+438, 0828+193 and 2104-242 we are 
confident that the seeing disk filled the slit. For 1809+407 
and 1931+480, it is reassuring that our results are in good 
agreement with Egami et al. (I2003II and Carson et al. (I2001II 
respectively. 

For the rest of the sample at least part ot the data where 
obtained with seeing size smaller than the slit width. 

We have estimated the maximum uncertainties on both 
the FWHM and Va values assuming the worse case scenario, 
i.e., the best seeing measured during the observations for 
each object and a clumpy morphology. The uncertainties 
estimated in this way are therefore upper limits. The ins- 
trumental profile for a spatially unresolved source would be 

seeinn size (arcsec) tti i tt-» • i • j i t 

. t 'j^i. r r X li , where IP is measured usmg the sky 

slit wzath {arcsecj ' o ^ 

lines (Table 1). The range of FWHM we obtain across the 
halos are: <550-790 km s"^ for 0211-122; 540-700 km s"^ 
for 0943-242; <710-960 km s"^ for 1410-001; 680-940 km 
s-i for 1558-003; and 700-810 for 2105+236. It could also 
be that some of the emission lines which we have found to 
be unresolved turn out to be resolved. 

We estimate that the maximum velocity shift (for Hell) 
which could be artificially introduced due to the spatial dis- 
placement between unresolved clumps in the direction per- 
pendicular to the slit would be ~310-365 km s~^ (depen- 
ding on the object) for 0211-122, 0943-242 and 1558-003 
and ~600 km s'^ for 1410-001 and 2105+236. 

Since all the estimations above have been done for the 
worse case scenario, we are confident that our main conclu- 
sion relative to the existence of quiescent gas in all objects 
are not affected. Our assumption that the objects filled the 
slit could have some non- negligible impact on some of the 



calculations in §4.4. This will be mentioned when appro- 
priate. 



4 DISCUSSION 

In this section we will investigate the origin of the velocity 
patterns, the ionizing mechanism, and the nature of the LS- 
BHs as well as the cosmological implications of our results. 

4.1 Nature of the kinematic components 

There is evidence that the kinematic and ionization proper- 
ties of high redshift radio galaxies can be determined in va- 
rying degrees by shocks driven during interactions between 
the radio structures and the ambient gas as well as pho- 
toionization by the central quasar (e.g. Best, Rottgering & 
Longair 2000). This is also the case of some intermediate 
redshift radio galaxies. As an example, Villar-Martfn et al. 
( 1999 ) explained the kinematics and optical emission line 
spectrum of the radio galaxy PKS2250-41 {z =0.3) in terms 
of two kinematic components one of which is emission from 
kinematically perturbed shocked gas and the other is non 
perturbed gas photoionized by the continuum from the ac- 
tive nucleus or the continuum generated by the shocked hot 
gas. 

We have found that all objects in our sample are asso- 
ciated with giant gaseous reservoirs with quiescent kinema- 
tics that sometimes extend beyond the radio structures. In 
addition, we have found highly perturbed gas in 8 out of 10 
objects in the sample, usually located inside the radio struc- 
tures. We propose that these are the emission from ambient 
non shocked gas and the emission from shocked gas respec- 
tively. 

There are other examples in the literature that 
support our interpretation. Spectroscopic observations of 
4C41.17 {z =3.8, van Breugel et al. 2001) and 1243+036 
l |van Ojik et al. 1996| l show that the gas kinematics is highly 
perturbed inside the radio structures and along the radio 
axis. Beyond the radio hot spot the velocity shifts and 
FWHM of the lines decrease abruptly. See also Maxfield et 
al. (.2002J . 

The discovery of a quiescent reservoir of gas in all the 
objects in our sample suggests that it might be a common 
ingredient of high redshift radio galaxies. 

4.2 Properties of the quiescent LSBHs 

4.2.1 Sizes 

Baum et al. (1988) and Baum & Heckman (1989) studied 
the emission line properties of a sample of 43 low redshift 
{z <0.6) powerful radio galaxies. In ~85% of these sources 
the line emission was spatially resolved. They found that 
the median total extent of the emission line nebulae is 10 
kpc, with only one object (3C227) with size (96 kpc) larger 
than the minimum total extension of the quiescent halos in 
our sample (~60 kpc; the maximum total extension is ~140 
kpc). The LSBHs are giant structures. 

There are a few examples in the literature of HzRG with 
Lya emission on the 100 kpc scale (McCarthy et al. 1990a, 
Francis et al. 2001, Pentericci et al. 1998, Reuland et al. 
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2003b). Our work shows that all the objects in our sample 
are associated with giant halos which extend far beyond the 
continuum structures revealed by the HST images and the 
high surface brightness emission line regions. 

4-2.2 Luminosities 

The Ha luminosities of the quiescent halos within the slit 
expected from the measured Lya luminosities (assuming 
case B recombination LyQf/HQ=8.8) are in the range ^0.1- 
2.8x10*^ erg s~^ (these are lower limits because Lya absorp- 
tion has been neglected). Baum & Heckman (1989) found 
that the median emission line luminosity in Hq:-|-[NII] in 
their low redshift sample is 3x10*^ erg s~^. Since the halos 
are likely to occupy a much larger spatial region outside the 
slit, it is clear that the total luminosities are probably much 
higher than line luminosities in low redshift radio galaxies. 

We have estimated the total Lya luminosity expected 
from the halos, assuming that the gas is photoionized by the 
quasar continuum (see §4.3.1) We emphasize the approxi- 
mate nature of these calculations due to the strong assump- 
tions relative to the geometry of the ionized gas, the neglec- 
tion of Lya absorption and the constant value of the Lya 
luminosity per unit volume across the halo. A biconical geo- 
metry with cone opening angle 90° is naturally expected in 
this scenario (Barthel 1989). We have assumed that the two 
cones have convex base and heights given by the extension 
of the LSBHs at both sides of the continuum centroid. The 
Lya luminosities (Table 2) and the volumes of the LSBHs 
within the slit and within the ionization cones were used to 
estimate an average luminosity per unit volume and the to- 
tal luminosity expected from the two cones. We obtain L^£ya 
in the range I.Txlff*^ - 1.2x10'*^ erg s~^. For comparison, 
Reuland et al. (2003b) measured total Lya luminosities of 
the giant nebulae in 3 radio galaxies at z ^3.4-3.8 in the 
range 5.9x 10*'*-1.4x 10*^ erg s~^. These values include the 
emission from the high surface brightness perturbed regions. 
This suggests that the quiescent gas emits a large fraction 
of the total line luminosity in HzRG. 

The implied Ha luminosities are > 8.7x 10*'^-1.4x lO*''' 
erg s~^, i.e., about > 100-1000 times higher than those of 
low redshift powerful radio galaxies (Baum et al. 1989). 

4-2.3 Kinematics 

The maximum velocity shift measured across the total ex- 
tension of the quiescent halos in our sample is <600 km s~^ 
and the maximum FWHM, except for 1410-001 and 1558- 
003, is ~700 km s"^ (Table 2). In those two radio galaxies, 
the maximum FWHM is ~800-850 km s^^ . It is important to 
note that spectral decomposition of Hell was not attempted 
in these two objects, but in both cases the presence of a 
narrower component is expected (see §3.1). 

Baum et al. (1990) studied the kinematic properties of 
the extended gas in a sample of 19 low redshift powerful 
radio galaxies. Emission line FWHM and maximum velocity 
shifts across the extended nebulae are in general < 500 km 
s~^. At low redshift, although rare, there are some objects 
where velocity shifts as high as ~800 km s"'^ have been 
measured and which present regular motions consistent with 
rotation (e.g. Tadhunter, Fosbury & Quinn 1989). 



Although the lines in the HzRG halos are somewhat 
broader the FWHM values are not extreme and do not im- 
ply kinematic perturbation. FWHM larger than 600 km s~^ 
are very rare in the extended nebulae (Baum et al. 1990) 
of low redshift radio galaxies, but some examples with no 
signs of kinematic perturbation exist: for instance, PKS 
13454-125, which was classified by Baum, Heckman & van 
Breugel (1992) as a 'rotator'. 

The broader lines in our sample might also be a conse- 
quence of the larger linear size covered by the apertures used 
in our kinematic analysis (>10 kpc, compared with <1 kpc 
in Baum, Heckman & van Breugel (1990) study). Velocity 
shifts >200 km s~^ across less than 10 kpc are common in 
low redshift radio galaxy nebulae. Such shifts will broaden 
the line profiles when integrating over larger (linear scale) 
apertures. As a reference for comparison, the typical line-of- 
sight stellar velocity dispersion of a giant elliptical galaxy 
is uios ~300 km s~^, corresponding to a characteristic line 
width (FWHM) of ~700 km s~^ (e.g. Heckman et al. 1991a). 

Therefore, the velocity shifts across the quiescent ha- 
los are consistent with measurements in low redshift radio 
galaxies. The lines are often somewhat broader, but the 
FWHM values are not extreme and consistent with quiet 
(non perturbed) kinematics. 



4-2-4 Densities 

Electron densities Ue for the quiescent halos were estimated 
using the relation Lhya = 4 x lO"^** nl f V erg s~^, where 
/ is the filling factor, and V the volume of ionized gas inside 
the slit (McCarthy et al. 1990a) The same geometry as in 
§4.2.2 was assumed: V is the slice of the convex cones in- 
cluded in the slit. Large uncertainties affect the calculations 
due to the assumptions on the geometry, the neglection of 
Lya absorption effects, the uncertainty on the / values and 
the assumption that the density is constant across the halo. 
/ values have been estimated to be ~10~^ ~10~* for typical 
emission line regions of HzRG dMcCarthy igMj l However, 
/ ~ values in the range 10~^-10~* were obtained in stud- 
ies of Lya nebulosities associated with high redshift quasars 
(Heckman et al. 1991b). The advantage in this case is that 
one knows the properties of the ionizing source with some 
confidence. 

We obtain densities in the range ~17-150 cm^'^ for 
/ =10"^ (a factor of 10 higher if /=10"''). In spite of 
all the sources of uncertainty, the calculated densities are 
similar to estimations by other authors using other tech- 
niques (e.g. Heckman et al. 1989, Heckman et al. 1991b) for 
gaseous nebulae associated with high redshift radio galaxies 
and quasars. As an example, McCarthy et al. ( 1990a I cal- 
culated an upper limit for the density of the ionized gas in 
3C294 of ~150 cm~^ based on estimations of the pressure 
in the diffuse radio plasma. 



* / values in low redshift radio galaxies (where the density can be 
constrained using the [SII] doublet) are in the range 
(Heckman et al. 1984, van Breugel et al. 1985) 
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4.2.5 Masses 

The masses are calculated as M = Ue V f nip, where nip 
is the proton mass and V is the volume within the ioniza- 
tion cones. For the densities derived above and assuming the 
same biconical geometry and / =10~^, the mass of ionized 
gas in the quiescent halos is in the range 10^-1° M0 (a factor 
of 10 lower for /=10"'^). For comparison, McCarthy (1993) 
quoted typical masses of ionized gas in HzRG of ~10^ Mq. 



4.3 LSBH energetics 

We investigate in this section the excitation mechanism of 
the giant LSBHs: photoionization by the quasar continuum, 
photoionization by the continuum generated by the hot 
shocked gas, photoionization by stars, cooling radiation from 
a cooling flow nebula. 



4.3.1 Quasar 

We estimate in this section the total quasar ionizing lumi- 
nosities Q\on implied by the Lya luminosities of the quies- 
cent halos. The total Lya luminosities expected from the 
ionization cones were calculated in §4.2.2. We have esti- 
mated Q\%n assuming that in case B recombination 68% of 
hydrogen ionizations yield a Lya photon (Osterbrock 1989). 
We have also assumed a covering factor C =0.15 of the quies- 
cent halos as seen by the nucleus (Heckman et al. 1991a). 
We obtain Q\l\ in the range 3.1x10^^ -7.0x10^*^ photons 
s-i. 

If the ionizing continuum follows a power law of index 
a=-1.0 ( Vill ar-Martm, Tadhunter fc Clark 1997| l, the total 
ionizing luminosities are in the range erg s 

In spite of the large uncertainties involved in the calcula- 
tions, it is interesting to note that the values we have ob- 
tained are in the same range as radio loud quasars (McLure 
& Dunlop 2001) at comparable redshift ^ . 

We conclude that the halo luminosities and measured 
sizes can be naturally explained by quasar photoioniza- 
tion. This is further supported by the emission line spectra 
(strong CIV, Hell and in some cases NV emission) typical of 
active galaxies, and inconsistent with star forming objects. 



4.3.2 Shocks 

We have calculated the UV ionizing luminosity Lfl^°'^'' 



ge- 



nerated by shocks of different velocities (range 500-900 km 
s~^, Bicknell et al. (12000')) which can explain the measured 
Lya luminosities LLya of the LSBHs ® . Combining equations 
(4.1) and (4.4) in Dopita & Sutherland (1996) and assuming 
case B recombination value for Lya/H/3=26.5: 



Lt^:^' = 4.2 X ( 



Vs 



X L 



Lya erg s 



100 km s-i' 

The total luminosity radiated by the shock is 

^ High redshift quasars {z =2-3.4) with luminosities as high as 
10** erg s~^ have been reported (Kaspi et al. 2003) 
^ The quiet kinematics of the halos suggests that this gas has not 
been shocked, therefore, we reject coUisional ionization. 



~2xLf^°'='= (Dopita & Sutherland 1996) and a large frac- 
tion should be emitted within the optical filters. Using the 
Lya luminosities measured within the slit for the objects 
in our sample the resulting Ll]^°'^^ values are in the range 
^several x IQ'**"*^ erg s~^ for 500 km s~^ shock velocity 
(higher values are obtained for higher velocities). 

There is, however, no evidence for such source of con- 
tinuum in the extended regions of HzRG. The high velocity 
gas gives us an idea of the spatial extent of the shocked gas. 
According to this, the continuum should have surface bright- 
ness in the range -several x 10"^'*-10"^'^ erg s ^ cm ^ 
arcsec"^. This is a gross lower limit, since we have only 
considered the gas within the slit. Typical optical surface 
brightness of HzRG measured from WFPC2/HST images 
are however < few x 10"^'' erg s~^ cm~^ arcsec"^. 

Another argument against shock processes is the ex- 
treme mass flow rates (M) through the shock required to 
explain the Lya luminosity (Villar-Martm et al. 1999). The 
required values are in the range 10^"" M© yr"^ for 14=500 
km s"'^ (a factor 2.1 smaller for Vs=900 km s~^). If the age 
of the radio source is ~10'^ yr and it is entraining mate- 
rial during its pass through the gaseous environment, this 
would imply that the total amount of material 'consumed' 
by the shock is in the range 10^°"" M©. These masses are 
larger than the total mass of ionized gas in high redshift 
radio galaxies (e.g. McCarthy et al. 1990a, van Ojik et al. 
1997, Villar-Martm et al. 2002). The amount of 'consumed' 
gas is likely to be higher, since we are only considering the 
one required to produce the luminosity measured inside the 
slit. 

It is therefore unlikely that the excitation of the quies- 
cent halos is due to shock related processes. 



4.3.3 Cooling radiation 

Steidel et al. (2000) proposed that the Lya giant ( > 100 
kpc) nebulae found around Ly break galaxies at 2 ~3 can be 
powered by radiative cooling of gas. Could this be the origin 
of the emission from the quiescent LSBHs in our sample? 

Models of high redshift cooling flows for galactic sys- 
tems in the process of formation lIFardal et al. 200T1 predict 
Lya luminosities due to cooling radiation as high as few x 
10^* erg s~^. According to these models the range of mea- 
sured Lya luminosities from the quiescent halos (Table 2) 
imply star forming rates (SFR) in the range ~50-700 M©. 
These are gross lower limits, since we have used only the 
luminosity within the slit and neglected Lya absorption ef- 
fects. 

Using the luminosity (within the slit) at 1500 A Vernet 
et al. (I2OOIII estimated a range in SFR of 2-60 Mq yr"^ 
for the objects in our sample, not taking into account dust 
reddening. Although the estimations are affected by large 
uncertainties, it is interesting to compare with the model 
predictions by Fardal et al (2001). If the ratio of un-obscured 
to obscured SFR is similar to that of Lyman break galaxies 
at z <3 (i.e. ~25, Sawicki & Yee 1998), the SFR in our 
objects could be in the range 50-1500 M© yr~^. If this is the 
case and according to Fardal et al. models, it is possible that 
in some cases the measured Lya luminosities are consistent 
with cooling radiation. 

A problem for these models is that the line emission is 
expected to be concentrated in a very small region of ra- 
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dius < 10 kpc or less for the most luminous halos (>few x 
lO'*'^ erg s^^). To explain the large sizes of the Lya nebulae 
around Ly break galaxies discovered by Steidel et al. (2000) , 
the authors propose that the Lya photons from gravitational 
cooling are scattered by neutral hydrogen in the outer neu- 
tral regions. However, we reject this interpretation since we 
have detected Hell and CIV in the outer parts of the quies- 
cent halos and therefore, this gas is ionized. 

Similar cooling flow models by Haiman, Spaans & 
Quataert I2000II predict Lya luminosities in the range 
-|^g43-44 gj.g g-i cooling radiation within measured 

angular sizes in the range 11 to 15 arcsec for z =2 to 3, 
consistent with the LLya values in Table 2. One of the im- 
portant differences with Fardal et al. (I2UU11 models is the 
assumption that Lya is produced out to the virial radius, 
while Fardal et al. (2001) predict a much more concentrated 
Lya emission. A problem for Haiman, Spaans & Quataert 
1200011 models is that the expected surface brightness for the 
halos is ^few x 10~^^ erg s~^ cm~'^ arcsec"'^, a factor of at 
least 100 lower than we observe. 

In summary, it is unlikely that the line emission from 
the giant quiescent halos is due to cooling radiation. Some 
high redshift cooling flow models can explain the measured 
Lya luminosities within the slit, however the models are in- 
consistent with the measured surface brightness or the an- 
gular sizes of the halos. In addition, a cooling flow nebula 
is likely to occupy a much larger volume than that within 
the slit and in such case, the high expected Lya luminosities 
would be also difficult to reconcile with cooling flow models. 

4.34 Stars 

Continuum has been detected from some of the LSBHs, 
which could be due to stars. We cannot confirm this, how- 
ever. Scattered light or nebular continuum might be domi- 
nant (Vernet et al. 2001). A strong argument against stellar 
photoionization is the detection of high ionization lines such 
as CIV, NV and Hell, which require a much harder ioni- 
zing continuum. The LSBH emission line spectra are typical 
of high redshift active galaxies and very different from high 
redshift star forming objects. 

Therefore, we conclude that the ionization of the quies- 
cent giant halos along the radio axis is dominated by the 
quasar continuum. We do not reject the possibility that 
other mechanisms excite the gas in other regions across the 
LSBHs. The giant Lya nebulae revealed by narrow band 
images obtained by Reuland et al. (2003b) show Lya emis- 
sion in regions which are also outside any putative ioniza- 
tion cones. Unless Lya is scattered (Villar-Martm, Binette 
& Fosbury, 1996) another mechanism must be responsible 
for the excitation of this gas. 

4.4 Origin of the halos and dynamical masses of 
HzRG 

4.4.1 Are the halos settled in a stable configuration ? 
Rotating disks 

V096 and VM02 proposed that the giant quiescent halos 
could be settled in rotating disks. Both observations and the- 
ory (IBarnes 200211 indicate that much of the interstellar gas 



in merging galaxies may settle into gaseous disks which may 
extend to several times the remnant half-light radii. Rota- 
tion of the extended emission line gas in some powerful radio 
galaxies at low redshifts suggests the presence of young disks 
of gas acquired in a recent interaction or merger with a gas 
rich galaxy ( |Baum, Heckman fc van Breugel. 1992^ . Large 
scale (up to 100 kpc) HI disk like structures have been found 
in several low redshift radio galaxies (Morg anti et al. 2002} 
and elliptical galaxies (O osterloo et al. 20Q 2|l. 

We have investigated whether there are signs of rotation 
in the LSBHs of our sample of HzRG. 

• 3 objects (0211-122, 1558-003 and 1931+480) show ve- 
locity curves inconsistent with rotation. 

• 4 objects (0943-242, 1410-001, 2104-242, 2105-f236) 
show velocity fields consistent with rotation curves with the 
slit located along or close to the line of nodes. For 1410-001, 
the velocity shifts at the positive side of the spatial zero 
might map half of the rotation curve (Fig. 4, bottom) 

• 1 object (1809-1-407) shows a velocity field consistent 
with a rotation curve with the typical V shape expected 
when the slit is located perpendicular to the line of nodes 
and with some impact parameter relative to the rotation 
center. This is also the case of 0828-f 193 (VM02) 

• for 0731-1-438 it is not possible to say, since the quiescent 
gas is detected only in two apertures. 

Therefore, 6 (including 0828-1-193) out of 10 objects in 
the sample have velocity fields consistent with rotation. The 
rotation center (probably located at the position of the ac- 
tive nucleus) is clearly shifted from the spatial continuum 
centroid in some objects. This is expected, since the opti- 
cal continuum is not direct AGN light, but nebular, stellar 
and/or scattered quasar light. 

When the rotation curve is measured along the line of 
nodes, the dynamical mass can be calculated as Mdyn ~ 
. , where R is the radius of the disc, V is half the 
amplitude of the rotation curve and i is the inclination an- 
gle of the disk with respect to the plane of the sky. We 
have calculated dynamical masses for 0943-242, 2104-242 
and 2105-1-236. R is the maximum extension of the LSBH 
from the rotation center. This was determined for 0943-242 
from the rotation curve (Fig 9, bottom). We have assumed 
that this position is marked by the radio core in 2104-242 
and 2105+236. This is a good assumption for 2105+236, 
since the radio core was aligned with the vertex of the ioni- 
zation cones (§2.2), where the active nucleus is expected to 
be. For 2104-242 the assumption might be wrong, since the 
NIR continuum centroid (which was aligned with the radio 
core) might be spatially shifted from the AGN. The results, 
however, are not strongly affected when considering a rea- 
sonable range of possible positions for the rotation center. 

We obtain Mdyn x sin^i ~2.8, 0.3 and 2.9x10^^ M© 
for 0943-242, 2104-242 and 2105+236 respectively^. For 
comparison, v096 estimated Mdyn x sin'^i ~10^^ for 
1243+036. Typical masses of cD galaxies at low redshift are 
~10^"^ M0 or more. If the radio galaxies above have similar 



^ If the gas across 0943-242 and 2105+236 did not fill the slit, the 
velocity curves might partially reflect the gas distribution within 
the slit rather than rotation (see §3.2). 
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masses, this implies disk inclination angles i < 30° relative 
to the plane of the sky. 

Alternative possibilities that would produce shallow ro- 
tation curves are 1) locating the slit across the rotation cen- 
ter and at certain angle (<90 degr) relative to the line of 
nodes 2) the objects are less massive than cD giant ellipti- 
cals at low redshift. This result was also obtained by Tad- 
hunter, Fosbury & Quinn ( iggg") from the kinematic study 
of the extended gas in low redshift powerful radio galaxies. 

Companion satellites 

As VM02 proposed, the halos might consist of a 
number of satellite companions. There is evidence that 
powerful radio galaxies at low redshift lie in rich environ- 
ments, with numerous companions at small distances from 
the main galaxy (within the inner 100 kpc). Some exam- 
ples are SClTi^^JHecknmi^^a^ Breugel fc Miley 1984| l, 
PKS2250-40 IITadhunter et al. 1994t aiid Coma A 
UTadhunter et al. 20001 . HzRG are known to lie in 
rich environments (Pentericci et al. 2002a, Venemans et al. 
2002) and numerous star forming companions have been 
found in the field of at least one high redshift radio galaxy 
JPentericci et al. 2002all . 

However, the continuous detection of the quiescent ha- 
los across tens of kpc and the AGN type spectrum suggests 
that they are gaseous reservoirs, rather than companion 
(proto)galaxies. Although we could think of a scenario where 
the central (proto)galaxy is surrounded by a high number of 
active galaxies l|Pentericci et al. 2002b1 or that the central 
quasar ionizes the gas in the companion galaxies, a more nat- 
ural explanation is that the halos are gaseous reservoirs. The 
diffuse and chaotic (filaments, plumes, cones) morphologies 
of the giant Lya nebulae revealed by recent narrow band im- 
ages of high redshift radio galaxies (Rculand et al. 2003bll 
further support this interpretation. Companion satellites 
might be embedded in the diffuse halos but the line emis- 
sion is likely to be dominated by the halos. Reuland et al. 
Il2003bl l find ' gaps' in the giant Lya nebulae which coincide 
with companion satellites revealed by K band images. 

Gaseous spherical envelopes 

Rather than settled in a disk, the quiescent halos could 
be gaseous envelopes supported against gravitation by the 
velocity dispersion within the halos (rather than rotation). 
The dynamical mass enclosed within a radius R can be cal- 

5 R 

culated as Mdyn ~ — cT^' where Vr is in this case the ra- 
dial velocity dispersion of the clouds llCarroll fc Ostlie 199911 
within the halo (given by the a of the emission lines) and R 
is Rmax (Table 2) . The estimated masses are shown in Table 
3. 

Therefore, Mdyn ^ several x 10^^ M© for most objects* . 
If the halos are virialized systems of clouds, these results sug- 
gest that most radio galaxies in this sample are less massive 
than nearby cD giant elliptical galaxies. 



* The calculations for 0211-122, 0943-242, 1410-001, 1558-003 
and 2105+236 could be affected by seeing effects (see §3.2). Con- 
sidering the worse case scenario the masses would change to 8.1, 
7.1, 12.1, 13.0, 8.5 Xl0l2, Mq respectively, using the maximum 
FWHM measured across the quiescent halos (Table 2) 



dyn 

hin RMax 



Name 


within RMax 


within RMax 


1809+407 


2.1 


4.1 


0211-122 


<1.2 


6.3 


1931+480 


2.3 


5.8 


1410-001 


<2.9 


8.3 


0731+438 


1.8 




2105+236 


3.2 


5.4 


2104-242 


<1.0 




1558-003 


4.4 


10.0 


0828+193 


<1.3 




0943-242 


2.3 


4.4 


1243+036 


1.0 





Table 3. Dynamical masses of high redshift radio galaxies as- 
suming that the gas consists of a virialized system of clouds. The 
second and third columns give the masses using the minimum 
and maximum FWHM of Hell across the halos (Table 2). Masses 
are given in units of lO^'^ Mq. Upper limits are given for those 
objects for which FWHMmi„ are upper limits. 



4.4-2 Inflows or Outflows 

The possibility that the quiescent halos are cooling flow 
nebulae was discussed in §4.2. Although we concluded that 
the line emission is likely to be powered by the quasar con- 
tinuum, this does not exclude the possibility that the halos 
are cooling flow nebulae (but with the quasar dominating 
the emission line processes). This will be discussed in more 
detail in §4.5. 

An alternative possibility is that the halo has been de- 
posited by galactic superwinds. Winds are needed to explain 
the chemical enrichment of the quiescent halos at tens of kpc 
from the nuclear region (the strength of the NV line rela- 
tive to CIV and Hell in 0211-122 and 0943-242 suggests at 
least solar metallicities, Vernet et al. 2001; Humphrey et al. 
2004 in prep.). We have calculated some basic parameters 
(following Heckman et al. 1990) characterizing a superwind 
which explains the kinematic properties and sizes of the 
quiescent giant halos: (1) energy injection rate dE/dt; (2) 
dynamical time tdyn', (3) infrared luminosity expected from 
the starburst Ljr ; (4) required star forming rate SFR; 
(5) supernova rate SNR; (6) mass injection rate dM/dt. 
The results for all the objects are shown in Table 4. We 
have assumed in our calculations a pre-shock density of 1 
cm~'^ (Heckman et al. 1990). Another parameter needed in 
the calculations is v, which is half the separation in velo- 
city between the two components of the double peak profile 
produced by the expanding bubble. The double peak is not 
resolved in our spectra and we can only estimate an upper 
limit for v. We estimate that two unresolved emission lines 
separated by a given velocity Sv would produce a line of 
FWHM>25v = 4v. Therefore v < EMHM^ where we have 
assumed FWHM = FWHM™i„ in Table 2. The values for 
dE/dt, LiR, SFR, SNR and dM/dt are likely to be high 
upper limits (lower limits for tdyn) due to the strong depen- 
dence of these parameters with v (oc u^). Typical v values of 
low redshift starburst driven galactic superwinds are seve- 
ral hundred km s~^ at several kpc from the center, v values 
might be smaller for the LSBHs, since the putative bubble 
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(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


Name 


dE/dt 


^dyn 


log(LiH/L0) 


SFR 


SNR 


dM/dt 




X 10'^'^ erg 


X 10* yr 












< 


> 


< 


< 


< 


< 


1809+407 


0.9 


3.2 


13.1 


3450 


26 


530 


0211-122 


1.1 


6.1 


13.2 


4200 


32 


650 


1931+480 


1.2 


4.7 


13.2 


4550 


34 


690 


0731+438 


1.2 


2.3 


13.2 


4570 


35 


700 


2104-242 


0.4 


9.8 


12.7 


1370 


11 


710 


1558-003 


4.8 


4.6 


13.7 


17970 


138 


2760 


0828+193 


0.6 


9.2 


12.9 


2200 


17 


340 


0943-242 


1.0 


5.9 


13.2 


3830 


29 


590 


1243+036 


1.0 {*) 


3.8 (*) 


13.2 


3000 {*) 


29 


255 


IRAS 00812-7112 


0.6 




12.9 


2000 




300 



Table 4. Outflows as the origin of the LSBHs. (2) Energy injection rate (3) dynamical time; (4) infrared luminosity expected from the 
starburst; (5) expected star forming rate; (6) supernova rate; (7) mass injection rate. (*) Results from v096. > or < indicate lower or 
upper limits respectively. The values for IRAS 00812-7112 have been taken from Heckman et al. 1990. 



has expanded to much larger distances (tens of kpc) and is 
likely to have decelerated. 

The derived Lm upper limits are consistent with those 
measured for ultraluminous infrared galaxies (ULIRGs, 
Lin ^W^^Lq). Similar rest frame FIR luminosities have 
also been estimated for high redshift radio galaxies from 
submm observations (Hughes, Dunlop & Rawlings 1997). 

Since the calculated parameters are gross upper limits, 
we can only say that a scenario such that the halos have been 
deposited by galactic winds is not excluded by the data set. 



4.5 Cosmological implications 

The discovery of the giant quiescent halos in all objects in 
the sample we have studied suggests that they are a com- 
mon ingredient of high redshift radio galaxies. Although the 
nature of the halos - be they disks, spherical envelopes, cool- 
ing flow nebulae, gas deposited by winds - is not clear, they 
are giant gaseous reservoirs within which the radio galaxy is 
embedded. 

Our results suggest that the halos surround the object 
completely. The gas we see in emission is the fraction of the 
halo inside the ionization cones which has been ionized by 
the quasar continuum. The gas outside these cones is likely 
to have much lower ionization level. Evidence for extended 
reservoirs of neutral gas associated with distant radio gala- 
xies has been found by van Ojik et al. (1997; see also de 
Breuck et al. 2003). 

We proposed in VM02 that the LSBHs could be part of 
the initial gas reservoir from which the galaxy started to con- 
dense. According to galaxy formation models (e.g. Haiman, 
Spaans and Quataert 2000), dark matter halos collapse in 
a cooling flow manner. As a result the baryons settle into 
rotationally supported exponential disks embedded in the 
dark matter halo. Circular velocities of ~few hundred km 
s~^ (consistent with our observations) would be expected 
for halos with masses > several x 10^^ Mq. 

According to Haiman, Spaans and Quataert (2000), the 
expected line widths are of order ~30[(l+2)/6] A , which 
corresponds to ~915 km s~^ for Hell. The observed FWHM 
of the lines from the quiescent halos can be easily reconciled 



with the models taking into account that most disks are 
probably not seen edge on. 

According to the models, due to Lya cooling radiation, 
the halos would have surface brightness ~few x 10~^^ erg 
cm~^ s~^ arcsec"^ out to the virial radius (>15 arcsec). 
However, if a quasar switches on at some point and illumi- 
nates the halos, the expected surface brightness would rise 
to r^lO"^**"^^ erg cm"^ s"^ arcsec"^ IHaiman fc Rees 200lt 
within the Strogrem radius consistent with our observations. 
The line widths, velocity shifts, measured angular sizes, line 
luminosities and surface brightness of the giant quiescent 
halos are consistent with this scenario. 

The interpretation of a cooling flow and a rotating disk 
do not exclude each other, they could be two consecutive 
phases of the same process. What is clear in our case, is 
that the halos are not made of pristine material. They have 
been already enriched with heavy elements, which suggests 
an advanced state in the process. We have also found pre- 
liminary evidence for rotation in some of the halos (although 
it is not possible to discard other interpretations). 

If this galaxy formation scenario is valid, we expect that 
similar structures exist associated with normal and active 
galaxies at different redshifts. The halos are likely to be more 
easily detected around powerful active galaxies, since the ac- 
tive nucleus provides an intense supply of ionizing photons. 
In fact, giant quiescent (total extension >100 kpc) emission 
line halos have been found in some z<l radio galaxies (e.g. 
Stockton, Ridgway & Kellogg 1996, Solorzano-Iharrea, Tad- 
hunter & Axon 2001). 

The quiescent LSBHs could be related (progenitors?) 
to the giant HI disk like structures found around some low 
redshift radio galaxies (Morganti et al. 2002^ and elliptical 
galaxies (Oostcrloo ot al. 2002). An important difference is 
that the HI absorption line in these objects is usually much 
narrower ( < 100 km s~^) than the emission lines from the 
quiescent LSBHs. This suggests that the high redshift struc- 
tures are less settled. 

It is possible that similar giant gaseous structures have 
also been found around non active galaxies at different red- 
shifts. Steidel et al. (2000) discovered two giant (physical ex- 
tent >100 h~^) diffuse Lya emitters apparently associated 
with previously known Ly break galaxies at < 2; >=3.09. 
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At low redshift, observations of absorption line systems in 
the spectra of background quasars have provided evidence 
for large {R ~100 kpc) extended gaseous envelopes that 
surround galaxies of a wide range of luminosity and morpho- 
logical type (e.g. Chen, Lanzctta & Webb 2001; Lanzotta ct 
al. 1995). These galactic envelopes might have similar origin 
and nature as the giant quiescent halos we have discussed 
around high redshift radio galaxies. 

Therefore, we propose a scenario such that a cooling 
flow has been triggered as a result of the collapse of a dark 
matter halo. The baryons will end up settling in a rota- 
tionally supported disk. During this process, the quasar has 
switched on and illuminates and excites the gas. This could 
happen before or after the gas has settled into the disk. The 
emission from the gas becomes easily detectable in this way, 
since its surface brightness is much brighter than expected 
from pure Lya cooling radiation. During the process, intense 
star formation has already taken place and outflows associ- 
ated with supernovae explosions or the quasar activity have 
enriched the halos with heavy elements at tens of kpc from 
the nuclear region. 



ture (rotating disks, spherical envelopes, cooling flow nebu- 
lae, gas deposited by galactic winds) the radio galax:y seems 
to be embedded within the halos. 

We propose a scenario such that the quiescent halos 
arc part of dark matter halos that collapsed in the early 
phases of formation of what will become a giant elliptical 
galaxy. A cooling flow is generated. For the more evolved 
halos, the baryons have settled in a rotationally supported 
disk. In other objects, the cooling flow is in an earlier phase. 
At some point, a quasar has switched on and the halos have 
become observable thanks to strong line emission powered 
by the quasar continuum. Intense star formation has already 
been triggered and associated stellar winds or quasar out- 
flows have enriched the halos with heavy elements at tens 
of kpc from the nucleus. The measured sizes, emission line 
widths, velocity shifts, presence of heavy elements, surface 
brightness and emission line luminosities of the quiescent ha- 
los arc consistent with this scenario. 2D spectroscopy is now 
critical to map the kinematic, ionization and morphological 
properties of the halos in two spatial dimensions. 



5 SUMMARY AND CONCLUSIONS 

We have studied the kinematic properties of the extended 
gas in a sample of 10 high redshift {z ~2.5) radio galaxies 
observed with the Keck II and VLT telescopes. We have dis- 
covered in all 10 HzRG (including 0828-^193, VM02) giant 
halos of quiescent gas with FWHM(HeII)<850 km s~^ and 
velocity shifts across the nebulae <600 km s~^. This com- 
ponent often extends over about 50 kpc from the nuclear re- 
gion and sometimes beyond the radio structures. The emis- 
sion line spectra of the quiescent halos are typical of active 
galaxies. In addition, 8 out of 10 objects contain kinemati- 
cally perturbed gas usually with high surface brightness and 
located inside the radio structures. 

We propose that we have isolated the emission from gas 
whose kinematics has been perturbed by jet-induced shocks 
and the emission from ambient non shocked gas (the quies- 
cent halos). 

The quiescent halos often extend for more than 100 kpc 
and sometimes beyond the radio structures. Typical Lya 
surface brightness and luminosities (within the slit) are in 
the range ~fewxlO~'^'^ *° erg s~^ cm~^ arcsec"^ and 
^IQis-ii gj.g g-i Egtimated densities are in the range ^^17- 
150 cm~^. The quiescent halos are enriched with heavy ele- 
ments at tens of kpc from the active nucleus. The quasar con- 
tinuum is the dominant excitation mechanism of the quies- 
cent halos along the radio axis. The implied quasar lumi- 
nosities are in the range ~severalx 10^^-10''^ erg s~^, in the 
same range as radio loud quasars at comparable redshift. 
Cooling radiation, shock photoionization and stellar pho- 
toionization cannot explain at the same time all measured 
properties (line luminosities, surface brightness, sizes, line 
ratios) of the quiescent halos along the radio axis. We do 
not reject the possibility that other mechanisms excite the 
gas in other regions of the halos. 

The quiescent halos are likely to be giant gaseous reser- 
voirs (rather than companion galactic satellites). The de- 
tection in all objects suggests that they could be a common 
ingredient of high redshift radio galaxies. Whatever their na- 
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Figure 1. 1809+407: Top: VLA radio map spatially aligned with the 2D Keck spectra of the main emission lines. The apertures 
selected for the kinematic analysis are also shown. The spatial zero in all figures is the position of the continuum centroid measured 
on the Keck spectra. Notice the LSBH extending beyond the radio structures with quieter kinematics compared to the high surface 
brightness regions. Bottom left; Spatially integrated spectra for the selected apertures. The flux scales are arbitrary to show clearly the 

emission lines. Bottom right: Kinematic properties of the ionized gas. FWHM and velocity shift relative to the Hell emission at the 
continuum centroid arc shown for the two kinematic components revealed by the spectral fit to Lya and Hell. Circles: Lya; triangles: 
Hell. Open symbols: narrow component; solid symbols: broad component. Arrows represent upper limits. The vertical lines mark the 
outer edge of the radio source. The narrow component has FWHM in the range ~500-700 km s~^ and maximum velocity shift of 400 
km s~^ across the nebula (as measured for HeH). The results are in good agreement with those obtained for Lya. 
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Figure 2. 0211-122: Symbols and lines as in Fig. 1. A very extended LSBH is detected in the 2D spectra of all emission lines (top). 
In all figures results for CIV are only shown when they provide convincing evidence for the presence of low velocity gas, a signature of 
the quiescent halos. Both CIV (squares) and Hell are very narrow in the outer regions of the object (FWHM<400 km s-\ bottom right 
panels). Hell has FWHM in the range <400-700 km s"! and maximum velocity shift of 200 km s ^ across the nebula. Lya is heavily 
absorbed at some spatial positions and the kinematic results obtained with this line are not reliable, especially in the inner apertures. 
Spectral decomposition of the lines was not performed for this object. 




Figure 3. 1931+480: Symbols and lines as in Fig. 1. A LSBH with apparently more quiescent kinematics compared to the high surface 
brightness regions is detected in the 2D Lya spectrum (top) . The fits to the Hell line reveal (bottom right panels) a narrow component 
with FWHM in the range ~450-720 km s'^ and maximum velocity shift of 450 km s across the nebula. Similar results are obtained 
for Lya. 
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Figure 4. 1410-001: Symbols and lines as in Fig. 1. An extended LSBH is detected in the 2D line spectra (top panels) with apparently 
more quiescent kinematics than the high surface brightness regions. The fits to Hell, Lya and CIV (bottom right panels) confirm the 
existence of narrow line emitting gas across several apertures with FWHM(Hell) in the range <472-800 km s"! and maximum velocity 
shift ~500 km s-i. 
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Figure 5. B3 0731+438: Symbols and lines as in Fig. 1. The double Lya peak in the spectrum of ap. 1 (bottom left panels) is due to 
the contribution of the separate component clearly seen in the 2D spectrum (top). We have not considered this as part of the quiescent 
LSBH. A narrow component is found across the rest of the object, with FWHM~600 km and velocity shift<100 km s"'^ (bottom 
right panel). The results obtained with the three lines are in good agreement. 
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Figure 6. 2105+236: Symbols and lines as in Fig. 1. Contrary to other objects, there is no evidence for a quiescent halo in the 2D 
line spectra of this object. However, the kinematic analysis (bottom right panels) reveals the presence of a narrow component with 
FWHM(HeII)~500-650 km s"^ and maximum velocity shift ~570 km s'^. 
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Figure 7. 2104-242: The spectra for this object were obtained at higher resolution (see text). Particularly interesting is the very narrow 
(FWHM<265 km s~^) Lya emission detected in the outer regions of the object (bottom right panels), clearly seen extending to the 
NE in the 2D Lya spectrum. This is the gas wc have considered as the quiescent LSBH. Two kinematic components are found in ap. 3 
to 5. We have represented them with black and grey colours (rather than solid and open symbols) to highlight the fa<;t that no broad 
components are found in this object. 
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Figure 8. 1558-003: Symbols and lines as in Fig. 1. Top: Notice the LSBH seen in Lya extending beyond the radio structures with 
apparently quieter kinematics compared to the high surfa<;e brightness regions. The kinematic analysis (bottom right) shows FWHM(Hell) 
in the range 550-850 km s"! and maximum velocity shift 550 km s ^ . There is some preliminajry evidence for more quiescent gas in this 
object (FWHM<410 km s'^, see text). 
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Figure 9. 0943-242: Symbols and lines as in Fig. 1. Notice the very extended LSBH with apparently quiescent kinematics extending 
towards the west well beyond the radio structures and in all the emission lines (top panels). A narrow component with FWHM~400-600 
km s~^ is detected in all apertures with maximum velocity shift across the nebula of ~450 km (bottom right). The results for CIV 
are in good agreement. 
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Figure 10. Spatially integrated spectrum of the quiescent LSBH 
in 0943-242 (emission detected beyond the radio structures). 
Compared with other objects, NV is very strong relative to CIV 
and Hell in this case, suggesting high levels of chemical enrich- 
ment. 



